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Abstiact 


A  computer  program  Is  developed  to  predict  the  performance  of  a  thermoacoustic 
generator.  The  defining  relationships  are  based  on  equations  derived  using  short  stack 
and  boundaiy  layer  approximations.  The  engine  modeled  Is  the  five  Inch  engine 
currently  in  operation  at  Los  Alamos  National  Laboratory.  The  working  fluid  In  the 
model  Is  a  helium-argon  gas  mix,  chosen  to  provide  flexibility  In  placing  the 
traiisuucer  relative  to  the  thermoacoustic  stack  assembly  while  maintaining 
resonance  near  the  desired  operating  frequency  of  120  Hz.  The  transducer  model  is 
based  on  a  linear  alternator  design  proposed  for  use  with  the  L.  A.  N.  L.  engine  in  the 
future.  The  thermodynamic  properties  of  the  gas  mix  are  determined  using  vlrial 
equations  and  blnarv-  gas  mixture  relationships.  The  defining  relations  are  solved 
simultaneously  and  an  Iterative  process  is  used  to  optimize  the  gas  mix  The  variables 
for  which  the  program  then  solves  are  pressure  amplitude,  hot  heat  exchanger 
teinperalure,  and  angular  frequency  as  functions  of  input  heat  flux,  system  geometry, 
and  load  characteristics.  From  these,  the  power  output  and  total  system  efficiency  are 
determined. 
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Complex  pressure  eimpUtudes 

Real  part  of  complex  pressure  amplitudes 

Sound  speed 
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Isobarlc  specific  heat  capacity 

Isochorlc  specific  heat  capacity 

Resonator  diameter 
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Geometry  factor 

Toted  energy  flux 

Electric  current 

Imaginary'  part  of 

Wave  number 

Thermal  conductivity 

Plate  half  thickness  or  Inductance 
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Pressure 

Heat  flux  per  unit  area 
Heat  flux 
Resistance 
Real  part  of 
Spring  constant 
Entropy  flux 
Time 
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Power  per  unit  volume 
Work  flux 
Position 
Plate  half  gap 

Impedance 

Imaginary  part  of  complex  pressure  amplitude  A 

Imaginary  part  of  complex  pressure  amplitude  B  or  thermal 

expansion  coefficient 

Normalized  temperature  gradient 

Specific  heat  capacity  ratio 

Penetration  depth 
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V  Kinematic  viscosity 
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Power  or  perimeter 
Density 

P^randtl  number 
Phase  angle 
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I.  INTRODUCTION 


A.  OVERVIEW 

A  great  deal  of  work  has  been  done  In  the  development  of  practical 
thermoacoustic  refrigerators  which  use  gas  as  a  working  fluid.  Elxamples  include  the 
Thermoacoustic  Refrigerator  (Hofler,  1986)  and  the  Space  Thermoacoustic  Refrigerator 
(Susalla.  1988).  Work  on  thermoacoustic  prime  movers  which  use  a  gas  as  a  working 
fluid  consists  primarily  of  enclosing  the  thermoacoustic  engine  In  a  tube  with  two  rigid 
ends  and  measuring  the  performance  of  the  engine  under  various  pressure,  temperature 
and  geometry  conditions.  Relatively  little  work  has  been  done  In  developing  a 
Ihermoacoustic  generator. 

The  history  of  thermoacoustics  Is  long,  but  significant  advances  in  practical 
thetmoacoustlc  devices  are  relatively  recent.  As  early  as  1777,  experlnients  were 
performed  by  B>Ton  Higgins,  in  which  acoustic  oscillations  were  excited  in  a  large  pipe 
by  suitable  placement  ol  an  internal  hydrogen  flame  IPutnain  and  Derniis,  1956). 
Hlggens'  work  eventually  led  to  the  modem  science  of  pulse  combustion,  whose 
applications  have  Included  the  Gemran  V-1  rocket  used  in  World  War  II  and  the 
residential  pulse-combustion  tumace  introduced  by  Lennox.  Inc.  in  198?  lSwifl.i988). 

The  earliest  antecedent  of  the  thermoacoustic  prime  movers  in  use  today  is  the 
Sondhauss  tube.  Over  100  years  ago,  it  was  observed  that  a  hot  glass  bulb  attached  to  a 
cool  glass  tube  sometimes  emitted  sound  at  the  tip  of  the  tube.  Sondhauss 
quantitatively  Investigated  the  relationship  between  the  pitch  of  the  sound  emitted 
and  the  dimensions  of  the  tube.  Lord  Rayleigh  explained  the  Sondhauss  tube 
qualitatively  In  1896.  but  a  quantitative  theoretical  description  of  these  phenomena 
was  not  achieved  until  about  20  years  ago  when  Rott  and  others  published  a  series  of 
papers  (Swlit,  1988).  Roll's  iheoiy’  combines  basic  principles  from  physics. 
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thermodynamics,  and  acoustics  to  describe  the  effects  found  in  both  thermoacoustlc 
prime  movers  and  heat  pumps. 

One  of  the  most  significant  advances  In  experimental  themioacousUcs  was  the 
realization,  by  Carter,  White,  and  Steele  of  North  American  Aviation,  Inc.  in  1962,  that 
placing  structures  (such  as  a  stack  of  plates)  inside  a  Sondhauss  tube  could  Improve  the 
performance  of  the  tube  In  generating  sound  (Swlfl,  1988). 

In  recent  years,  a  number  of  people  Including  Wheatley,  Hofler,  Swift,  and 
Mlgllorl,  performed  a  series  of  experiments  at  Los  Alamos  National  Laboratory 
(L.A.N.L.)  In  New  Mexico.  The  purpose  of  the  experiments  was  to  investigate  the 
thermoacoustlc  effects  in  both  prime  movers  and  heat  pumps,  and  compare  these 
experimental  results  to  Rott's  theory  (Susalla,  1988). 

The  first  acoustic  heat  pump  built  at  L.A.N.L.  consisted  of  a  speaker  mounted  at 
one  end  of  a  closed  tube  which  contained  a  slack  of  fiberglass  plates  positioned  near  the 
closed  end.  They  found  that  a  temperature  difference  of  1(X)  degrees  centigrade  across 
the  stack  of  plates  was  easily  produced  In  about  a  minute  with  the  speaker  driving  at 
acoustic  resonance  fWheatly,  Swift,  and  Mlgllorl,  1986) 

Shortly  after  the  first  heat  pump  was  built  at  L.A.N.L,,  Hofler  built  a  device, 
which  has  come  to  be  (mown  as  the  Hofler  tube,  to  use  as  a  demonstration  for  his  Ph.D. 
candidacy  committee  at  the  University  of  California,  San  Diego.  The  device  {Fig.  1) 
produces  an  impressive  level  of  sound  from  the  open  end,  about  lOOdB  at  one  meter 
(Swift.  1988). 

Other  devices  built  and  studied  at  L.A.N.L.  include  Hofler  s  refrigerator,  the  heat- 
driven  acoustic  cooler  (Fig.  2),  a  liquid  sodium  acoustic  engine,  and  a  liquid  propylene 
engine  (Whealfy,  Swift,  and  Mlgllorl,  1986).  Currently  being  studied  at  L.A.N.L.  Is  a  five 
Inch  thermoacoustlc  prime  mover,  which  will  be  described  in  detail  later.  It  is  the 
prime  mover  upon  which  the  calculations  described  heiin  are  based. 
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B.  MOTIVATION 

The  results  of  these  early  ee^perincnls  on  themioacoustic  refngerators  at  LA.N.L. 
indicated  a  system  performance  that  was  far  less  than  expected  (Holler,  1986).  These 
results  led  Hofler  to  apply  Roti's  theory  to  the  expenmenial  systems,  solving  the 
resulting  equations  numencally.  He  also  built  the  thermoacoustic  refrigerator,  made 
accurate  m.easurements  of  its  thermodynamic  efficiency,  and  coinpaied  these  results 
to  those  predicted  by  theory.  The  work  which  followed  at  the  Naval  Postgraduate 
School  by  Susalla  and  others  was  done  to  unprove  the  overall  efficiency  of  the  system 
and  make  it  suitable  lor  space  applications  (Susalla,  1988. 

Swift  s  research  group  at  LA.N.L.  is  currently  persuing  a  similar  line  of 
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ure  2:  The  heal  dnven  acousuc  cooler.  High-lemperaiure  heat  causes  low 
'.peraturc  refrigeration  (Swift,  1988) 


investigation  for  a  thermoacoustic  generator.  A  five  Inch  diameter,  four  meter  long 
prime  mover  has  been  built  and  operated  In  a  rigid-rigid  resonator.  The  system  is 
operated  at  a  mean  pressure  of  twenty  atmospheres  with  hehum  as  the  working  fluid. 
Acoustic  pressures  on  the  order  of  ten  percent  of  the  mean  pressure  have  been  observed. 
Acoustic  ;jower  generated  is  on  the  order  of  one  kilowatt. 

The  next  phase  of  the  develop  nent  of  an  acoustic  generator,  using  the  five  inch 
engine.  Is  to  obtain  a  transducer  able  to  convert  the  acoustic  energy  to  electrical  energy. 
Two  contractors,  Methanlcal  Technology  Inc.  of  Latham,  New  York  and  Sun  Power  Inc. 
of  Athens,  Ohio  have  been  requested  to  provide  quotations  for  deliver^’  of  a  linear 
alternator  for  this  purpose.  The  cost  of  such  a  device  has  been  estimated  at  between 
$2F,000  (Sunpower)  and  $1 15,000  (MTI).  Because  of  the  high  cost  (and  current  lack  of 
funding)  It  was  decided  that  a  computer  program  should  be  developed  In  the  interim, 
which  would  predict  the  performance  of  the  combined  prime  mover-transducer  system. 
In  addition,  such  a  program,  once  experimentally  validated,  could  be  used  to  determine 
the  effect  of  changing  the  geometry  or  material  of  the  stack  assembly,  the  position  of 
the  stack  assembly  In  the  resonator,  or  the  characteristics  of  the  mechanical  or 
electrical  load  on  the  performance  of  the  system  as  a  whole.  Such  a  program  could 
result  in  a  great  deal  of  sa\'ings  In  material  and  man  hours  when  design  modifications 
are  contemplated. 
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n.  THEORY 


A.  THERMODYNAMICS 

In  designing  any  heat  engine,  compromises  are  made  between  thermodynamic 
Ideals,  and  engineering  and  marketing  realities.  In  designing  any  device,  a  practical 
balance  Is  sought  between  the  desire  for  high  efriciency,  power,  and  reliability  and  the 
necessity  of  achieving  low  cost,  small  size  and  maximum  simplicity.  Thermoacoustic 
engines  are  extremely  simple  devices,  as  will  be  seen,  and  promise  to  deliver  a 
substantial  fraction  of  the  Ideal  thermoacoustlc  efficiency. 

In  any  basic  thermod>'namics  course  we  learn  that  there  are  two  principle  types  of 

heal  engines:  heat  pumps  and  prime  movers.  A  heat  pump  removes  heat  Qc  from  a  low 
temperature  reservoir-  at  temperature  Tc.  receives  workW  from  the  surroundings,  and 
delivers  heat  Qh  to  a  high  temperature  sink  at  temperature  TTh. 

The  process  Is  reversed  In  a  prime  mover.  Heat  Qh  Is  received  from  a  high 
temperature  reser%-otr  at  temperature  Th.  work  W  Is  done  on  the  surroundings,  and  the 

remaining  heat  Qc  Is  dumped  to  a  low  temperature  sink  at  temperature  Tc- 

Tlie  first  and  second  laws  of  thermodynamics  place  an  upper  limit  on  the 
efficiency  of  any  heat  engine.  This  limit  Is  best  understood  by  applying  the  Carnot  cycle 
to  the  process.  The  Carnot  cycle  was  first  introduced  in  1824  by  N.  L.  Sadi  Carnot  (Sears 
and  Salinger,  1975).  The  Carnot  cycle  (Fig.  4)  consists  of  two  reversible  isothermal 
processes  and  two  reversible  adiabatic  processes.  The  direction  of  the  cycle  (l.e. 
clockwise  or  counter-clockwise)  determines  whether  the  engine  is  a  prime  mover  or  a 
heal  pump.  Since  the  work  In  this  thesis  deals  solely  with  prime  movers,  the  remainder 
of  the  discussion  In  this  section  will  focus  on  them. 
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From  the  first  law  of  thermod>'namics  we  know  that  energy  must  be  conserved  In 


the  system.  It  follows  that 

Qh-9c  =  w.  (1) 

The  second  law  requires  that  the  total  entropy  In  the  system  must  Increase  or  remain 
constant.  The  change  In  entropy  for  Isothermal  heat  transport  to/from  a  resorvolr  is 
given  by 


////// 


/////////. 


Figure  3:  Thermodymamic  cycles  of  a  heat  pump  and  prime  mover. 


AS  =  Q/T. 


(2) 


For  a  prime  mover  the  second  law  may  be  written  as 

ASC  -  ASh  =  Qc/Tc  -  Qh/TIi  >  0.  (3) 

The  efficiency  of  interest  in  a  prime  mover  (the  work  output  divided  by  the  heat  input)  Is 

n=  W/Qh.  CombiniTig  Eqs.  1  and  3  it  can  be  shown  that 


T1  =  W/Qh  <  (Th  ■  Tc)/Th  (4) 

The  temperature  ratio  on  the  right  side  of  the  above  equation  Is  knowm  as  Carnot's 


efficiency  (tie  =  (Th  -  Tc)/Th)  and  is  the  highest  efficiency  that  a  prime  mover  can 


achieve. 
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a  THERMOACOUSTIC  THEORY 

As  stated  In  Chapter  I,  thermoacoustlc  theory  was  developed  In  detail  by  Rott  and 
others  in  the  late  1960's  and  1970's.  Since  then,  a  great  deal  of  work  has  been  done  by 
Swift,  Hofler,  Wheatley,  and  others  in  applying  Rott's  theory  to  experimental  devices. 
The  work  contained  in  this  thesis  relies  heax'lly  on  the  relationships  developed  by 
applying  the  short  stack  and  boundary  layer  approximations  to  the  most  general  of 
thermoacoustic  derivations  (Swift.  1988). 

L  A  gualltatlve  Picture 

Consider  a  gas  parcel  In  thermal  contact  with  and  oscillating  along  a  plate  ai 
the  acoustic  frequency  (Fig.  5).  As  the  gas  parcel  moves  back  and  forth,  it  will  experience 
changes  in  pressure,  volume,  and  temperature.  The  temperature  change  is  due  to  both 
the  adiabatic  compression  and  expansion  of  the  gas  by  the  acoustic  pressure  and  by  the 
local  temperature  of  the  plate. 

Consider  first  the  gas  parcel  in  motion  along  a  plate  with  no  temperature 
gradient  (Fig.  5a).  The  gas  parcel  starts  at  position  1  with  a  temperature  T.  The  acoustic 
wave  moves  the  parcel  to  the  right  a  distance  dx  to  position  2  and  compresses  it 

adlabatlcally,  raising  Its  temperature  to  T*"*’.  The  temperature  of  the  gas  Is  now  higher 
than  the  local  temperature  of  the  plate  so  heal  Q  flows  from  the  gas  to  the  plate,  causing 

the  plate  to  warm  up  at  position  2  and  the  gas  to  cool  to  temperature  T^.  As  the  gas  cools 
Its  volume  decreases,  so  that  work  is  done  on  the  parcel  by  the  surrounding  gas.  The  gas 
parcel  now  moves  back  to  position  1,  undergoing  adiabatic  expansion.  It  arrives  at 

position  1  with  a  temperature  T.  which  Is  lower  than  the  local  plate  temperature,  The 
result  Is  a  heat  flow  from  the  plate  to  the  gas  parcel.  The  increase  in  temperature  of  the 
gas  parcel  causes  It  to  expand,  resulting  in  w’ork  being  done  by  the  parcel  on  the 
suiTounding  gas.  As  this  cycle  is  repeated,  the  constant  dumping  of  heal  at  position  2 
and  the  removal  of  heat  from  position  1  results  in  a  temperature  gradient  across  the 
plate  and  heat  transfer  to  the  right  along  the  plate. 
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Next,  consider  the  cycle  for  a  prime  mover  (Fig.  5b).  Assume  a  temperature 
gradient  (VT)  is  maintained  along  the  plate  by  means  of  heat  exchemgers  located  at 
either  end  of  the  plate,  with  the  hot  heat  exchanger  located  to  the  right.  Now  when  the 
gas  parcel  Is  compressed  adlabatlCcilly  and  moved  from  position  1  to  position  2,  heat 

will  flow  from  the  plate  to  the  gas  parcel  as  long  as  T  +  dx(VTl  ■>  'T'"'".  This  process  will 

result  in  an  expansion  of  the  gas  parcel  as  it  Is  heated  to  a  temperature  so  that 
work  Is  done  on  the  surrounding  gas.  When  the  gas  Is  then  adlabatlcally  expanded  and 

moved  back  to  position  1  It  will  be  at  temperature  Heat  must  then  flow  from  the  gas 
to  the  plate  and  work  will  be  done  on  the  parcel  by  the  surrounding  gas.  Heal  Is  thus 
pumped  from  the  hot  heat  exchanger  to  the  cold  heat  exchanger. 

The  actual  process  is  sinusoidal  rather  than  articulated  as  described  above. 
As  we  shall  see  later,  the  heat  transfer  lakes  place  only  within  a  thermal  penetration 

depth  of  the  surface  of  the  plate.  This  thermal  penetration  depth  (5ic)  is  the  distance  heat 
can  diffuse  through  the  gas  In  a  time  equal  to  the  reciprocal  of  the  angular  frequency  at 
which  the  motion  Is  taking  place.  Thermoacoustic  engines  may  be  thought  of  as  a  long 
train  of  adjacent  gas  parcels,  all  within  about  a  thermal  penetration  depth  of  the 
surface  of  the  plate,  that  draw  heat  from  the  plate  at  one  extreme  of  their  motion,  and 
deposit  heat  at  the  other  extreme  (WTieatlcy,  Swift,  and  Migllorl.  1986). 

2.  The  Single  Plate 

The  basic  principles  of  thermoacoustic  engines  are  most  clearly  Illustrated 
by  examining  the  simple  case  of  a  single,  small  solid  plate  aligned  parallel  to  the 
direction  of  vibration  of  a  standing  acoustic  wave.  The  presence  of  the  plate  In  the 
standing  wave  modifies  the  wave,  resulting  in  (1)  a  time-averaged  heat  flux  near  the 
surface  of  the  plate,  along  the  direction  of  the  acoustic  vibration  and  (2)  the  generation 
or  absorption  of  real  acoustic  power  near  the  surface  of  the  plate.  These  two  simple 
effects  are  the  basis  of  all  thermoacoustic  engine  phenomena  (Swift.  1988).  The 
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following  discussion  follows  Swift's  paper  closely.  The  only  significant  difference  i.s 
that  we  have  chosen  to  model  the  acoustic  pressure  as  a  cosine  function  so  that  the 
resulting  expression  for  fluid  particle  velocity  contains  a  sine  term. 
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In  the  absence  of  the  plate,  the  pressure  and  x-dlrection  velocity  of  the 
standing  wave  have  well  known  acoustic  properties.  Let  the  fluid  pressure  be  given  by 

p  =  pm  +  pie**^  (5) 

where  the  acoustic  pressure  is  given  by 

pi  =  Po  cos(kx)  5pi®(x).  (6) 


Figure  6:  Geometry  of  the  single  plate  thermoaccustic  example. 


The  fluid  particle  velocity  Is  found  by  applying  the  linear  invlscld  force  equation 
(Kinsler  and  others,  1982)  to  Elq.  5.  The  result  is 

u  =  uie“^.  (7) 

where 

ui  =  -l(Po/pmc)  sln(kx)  s-lui®(x).  (8) 

In  the  absence  of  the  plate,  the  sound  wave  Is  adiabatic  and  has  an 
oscillating  temperature  given  by 

Tl  =  fTniP/pmCplpl.  (9) 

For  an  Ideal  gas  we  know 

TmP/pmCp  =  (7-l)Tm/YPm.  (10) 

Combining  Eqs.  9  and  10  it  can  be  shown  that 

Tl/Tm  =  ((‘i^ll/Ylpi/pin  (11) 

for  ideal  gases.  Thus  the  fractional  temperature  oscillations  are  of  the  same  order  of 
magnitude  as  the  fractional  pressure  oscillations. 

The  Introduction  of  the  plate  Into  the  standing  wave  modifies  the 
unperturbed  temperature  oscillations  due  to  the  heal  flow  between  the  gas  and  the  plate. 
The  presence  of  the  plate  inodiiles  the  temperature  oscillations  In  both  magnitude  and 
phase,  for  fluid  within  about  a  thermal  penetration  depth  of  the  plate  surface.  The 
thermal  penetration  depth  In  the  fluid  is  given  by 

8k  =  V(2k7u).  (12) 

At  this  point  It  Is  prudent  to  make  several  assumptions  that  will  facilitate 
the  calculations  which  follow.  These  are: 

(1)  that  steady  state  exists, 

(2)  the  plate  Is  short  compared  to  the  acoustic  wavelength  and  that  It  Is  located  far 
from  both  the  pressure  and  velocity  nodes  so  that  both  pi  and  ui  can  be  considered 
uniform  over  the  plate. 
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(3)  the  Quid  Is  invlscid  so  that  ui  does  not  depend  on  y. 

(4)  the  plate  has  a  large  enough  heat  capacity  per  unit  area  so  that  it’s  temperature 
does  not  change  appreciably  at  the  acoustic  frequency, 

(5)  the  thermophysical  properties  of  the  plate  material  and  the  gas  are 
I"  dependent  of  temperature. 

(l,  the  plate  has  a  mean  temperature  gradient  In  the  x  direction  VTm  and  the  plates 

thermal  conductivity  In  the  x  direction  Is  negligible, 

(7)  the  fluids  thermal  conductivity  in  the  x  direction  Is  also  negligible, 

18)  the  mean  fluid  temperature  Is  Independent  of  y  and  Is  the  same  as  that  of  the 
plate. 


Calculation  of  the  oscillating  fluid  temperature  is  accomplished  using  only 
the  first  order  terms  In  the  general  equation  of  heat  transfer,  expressing  the  entropy  In 
terms  of  the  pressure  and  temperature  of  the  fluid  as 

SI  =  (cp/Tm)Tl  -  ((i/piri)pi.  (13) 

and  finally  solving  the  resulting  differential  equation  subject  to  the  appropriate 
boundary  conditions  (Swift.  1988.  p.  1151).  The  solution  Is 


Ti=|(Tmp/pmCp)pis  -  (VTm/a))u  1^1(1  -  e-U+l)y/6K).  (14) 


To  understand  the  terms  of  this  equation,  first  lake  the  limtt  as  6k'  «  y  where  the  fluid 
makes  negligible  thermal  contact  with  the  plate.  Under  this  condition 


Tl  ^(TmP/pmCp)pi®  -  (VTin/wlui®.  (15) 

The  first  term  is  due  to  the  adiabatic  compression  and  expansion  of  the  gas.  The  second 
term  Is  a  result  of  the  mean  temperature  gradient  in  the  fluid  and  the  fact  that  the  fluid 

oscillates  in  the  x  direction  with  displacement  amplitude  ui^/w.  Setting  Ti  =  0  and 
solving  Eq.  15  for  VTm.  we  see  that  there  Is  a  critical  mean  temperature  gradient 


VTcrll  =  TmPMpi^/pmCpJl®  (16) 

for  which  the  temperature  oscillations  at  a  point  are  zero.  This  critical  temperature 
gradient  is  Important  because  It  marks  the  boundary  between  the  heat  pump  and  the 
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prime  mover  regions  of  operation.  Further  manipulation  of  Eq.  16,  using  Ekjs.  6  and  8. 
and  the  generalization  of  Eq.  10  to  an  arbitrary  fluid  leads  to  the  useful  result 

VTcrit  =  l(Y-l)/TmPl(Tmk)  cot(kx)  (17) 

The  exponential  term  In  Eq.  14  is  complex.  It  approaches  zero  for  y  «  5tc 

where  the  plate  Imposes  the  boundary  condition  Ti  =  0.  When  y  =  5k  it  has  a  substantial 
imaginary  part.  This  phase  shift  In  the  oscillating  temperature  of  the  standing  wave 
due  to  the  presence  of  the  plate  leads  directly  to  the  time  averaged  heat  flow  in  the  x 
direction. 

a  Heat  Flux 

Since  wc  have  assumed  that  thermal  conductivity  in  both  the  plate  and 
the  gas  is  negligible  in  the  x  direction,  the  only  way  that  heat  can  be  transported  In  the  x 
direction  is  by  the  hydrodymamlc  transport  of  entropy,  thus  the  heat  flux  per  area  is 
given  by 

q2  =  Tmpm<siui>.  (18) 

Substituting  Eq.  13  for  si  and  simplifying  results  in  an  expression  for  heat  flow  as  a 
function  of  y.  The  total  heat  flux  Q2  along  the  plate  in  the  x  direction  can  then  be  found 
by  integrating  q2  over  the  y-z  plane.  The  resulting  expression  is 

92  =  -(l/4)n6KTmPpiHil®(r-l)  (19) 

where  F  =  VTm/^Tcnt  is  the  ratio  of  actual  temperature  gradient  to  critical  temperature 

gradient.  Tlie  term  TmP  is  nearly  1  for  ideal  gases  and  Uie  term  HSk  Is  the  area  which  is 
thermodynamically  active  in  a  plane  perpendicular  to  the  plate.  The  heat  flux  will  be 
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zero  If  the  plate  Is  positioned  at  either  a  pressure  or  a  velocity  node  of  the  standing  wave 


which  Is  readily  understood  from  the  product  pi®ui®.  Finally  we  see  that  If  VTm=  VTcrli 

then  r  =  1  and  no  heal  flux  occures.  If  VTm  <  VTcrlt  .  then  F-l  <  0  and  the  heat  flux  Is 

toward  the  velocity  node.  If  VTm  >  VTcrt'. .  then  F-l  >  0  and  the  heal  flux  Is  toward  the 
pressure  node.  If  suitable  heat  exchangers  are  positioned  at  the  ends  of  the  plate,  the 
heat  flux  carries  heat  from  one  heat  exchanger  to  the  other. 

b-  Work  Flux 

The  work  flux,  or  acoustic  power  Is  derived  In  a  similar  manner  from 
thermodjmamic  principles.  The  work  done  by  a  differential  volume  element  of  fluid  in 
going  from  a  volume  dx  dy  dz  to  a  volume  dx  dy  dz  +  dV  Is 

dw  =  p  dV'’  =  •  (p/p)  dx  dy  dz  dp.  (20) 

so  the  power  per  unit  volume  Is 

w  =  -  (p/p)  dp/dl  (21) 

Since,  for  osclllatoiy  motion,  we  are  interested  in  the  time  average  work  flux  we  use  Eq. 
5  for  p  and 

dp/dt  =  lojpi  +  ul(^pm/^x)  (22) 

to  find 


w2  =  -  (oypm)  <lpipi>  (23) 

for  the  average  acoustic  power  produced  (or  absorbed)  per  unit  volume.  /Vfter  expressing 
pi  In  terms  of  Ti  and  pi  and  substituting  into  Eq.  22  we  Anally  arrive  at  the  result 


w2  =  -  (1/2)  uPpi®  ImlTl).  (24) 

Just  as  w'ltl’  the  heat  flux  per  area,  only  the  Imaginary  part  of  the  temperature 
oscillations  contributes  to  the  acoustic  power  produced  per  unit  volume.  The  fluid 
within  about  a  thermal  penetration  depth  of  the  surface  of  the  plate  "breaths"  because  of 
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thermal  expansion  or  contraction,  with  the  correct  time  phase  with  respect  to  the 
oscillating  pressure  to  do  or  absorb  net  work. 

By  Integrating  W2  over  all  space,  using  Eq.  14  for  Ti ,  the  total  acoustic  power 
produced  Is  found  to  be 

W2  =  (l/4)n6KAx(Tmp2a)/pinCp)(piS)2(r-l).  (25) 

This  acoustic  power  Is  proportional  to  the  volume  of  the  fluid,  FISkAx,  that  Is  within  a 
thermal  penetration  depth  of  the  plate.  U  Is  also  quadratic  In  the  acoustic  amplitude 
and  vanishes  at  the  pressure  nodes  Finally  ll  is  proportional  to  the  same  temperature 
gradient  factor  T-l  as  appeared  In  the  expression  for  the  total  heat  flux.  Whether  this 
power  Increases  the  amplitude  of  the  standing  wave,  Is  radiated  away  to  Infinity,  or 
flows  through  an  acoustlc-to-electrlc  transducer  to  generate  electric  power  depends  on 
the  details  of  the  resonator. 

c  Efficiency 

The  efflclencj'  of  the  single  plate  engine,  used  as  a  prime  mover  Is  given 
by  T]  =  W/Qh.  For  a  short  plate,  Q2  is  essentially  constant  with  x.  so  9h  =  Qc  =  Q2. 
Substituting  for  W2  and  Q2  from  Eqs.  19  and  25  above  results  in  the  expression 

T)=  AxPupi^/pmCpuiS.  (26) 

Using  the  definition  of  VTcrlt,  this  equation  can  be  recast  Into  the  form 

q  =  ric/r.  (27) 

Tlius  vvc  see  that  as  r-+  1,  the  single  plate  engine  approaches  Carnot's  efficiency,  but 
this  Is  also  the  point  wliere  the  power  output  of  the  engine  approaches  zero. 

3.  The  Short  Engine 

.1,  the  prev’iour  section,  thermoacoustic  effects  were  examined  for  the  simple 
case  of  f  rgle  plate  in  an  acoustic  standing  wave.  In  order  to  examine  these  elTects 
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without  a  great  deal  of  distractions,  several  slmpllf>lng  assumptions  were  made.  The 

central  ihermoacoustlc  effects  that  were  revealed  are: 

(1)  The  thermal  boundary  condition  imposed  on  the  fluid  by  the  plate  causes  a 
phase  shift  of  the  oscillating  temperature  of  the  fluid  within  a  thermal 
penetration  depth  of  the  plate. 

[Zt  Heat  flows  hydrodynamlcally.  parallel  to  the  plate  In  the  fluid  about  a  thermal 
penetration  depth  from  the  plate. 

(3)  Acoustic  power  is  absorbed  or  produced  by  the  fluid  about  a  thermal 
penetration  depth  from  the  plate. 

(4)  The  direction  of  heat  flow,  and  whether  power  is  produced  or  absorbed,  depends 
on  the  magnitude  of  the  mean-temperature  gradient  relative  to  the  critical 
temperature  gradient. 

We  now  turn  our  attention  to  the  more  realistic  model  of  an  engine  built  with 
a  stack  of  parallel  plates.  The  two  sets  of  coordinate  axis  used  In  the  detailed 
calculations  (Swift.  1988.  pp.l  176-1 179)  which  are  summarized  below  are  Indicated  In 
Fig.  7.  We  shall  also  discard  many  of  our  previous  simplifying  assumptions  although  we 
shall  retain  the  assumptions  that  the  temperature  spanned  is  much  less  than  the 
absolute  temperature  and  that  the  length  of  the  plates  is  much  less  than  an  acoustic 
wavelength. 

Swift's  derivation  proceeds  as  follows.  Consider  a  parallel  plate  engine 
geometry’  w’lth  the  acoustic  vibration  along  the  x  axis  and  the  y  axis  perpendicular  to 
the  planes  of  the  parallel  plates.  The  equation  of  motion  for  the  fluid  is 
p(()v/dt  +  (vV)v) 

=  -Vp  +  +  (4+p/3)V(V’v),  (28) 

with  the  boundary  condition  v  =  0  at  the  surface  of  the  plates.  The  continuity  equation 
for  the  fluid  is 


3p/9l  +  V-(prv)  =  0. 

The  equations  for  heat  flow  In  the  fluid  and  solid  are 


(29) 


pTOs/dt+v*Vs)  =  V(KVD 

+  (terms  quadratic  In  velocity)  (30) 
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and 


psCs3Ts/dt  -  KsV'^Ts, 


(31) 


Stack  cf  plates 


Solid  y 


r  1  u  1 C 


Fluid 


Solid 


Flfjure  7;  Geometry  used  for  mulllplate  engine  calculations,  (a)  0\-er^'iew  of  short 
engine  In  resonator,  and  fb)  expanded  view  of  plates.  Each  plate  has  thickness  2L 
and  Is  separated  from  adjacent  plates  by  a  distance  2yo. 


p  =  Pm  +  piWe**^. 


(32) 


p  =  pmlx)  +  pilx.y)e‘“^  (33) 

V  =  lui(x.y)e‘“^+  Jvi(x.y)e‘‘^,  (34) 

T  =  Tm(x) +Ti(x.y)‘“‘.  (35) 

Ts  =  Tm(x) +Tsl(x.y)“^.  (36) 

s  =  Sm(x)  +  si(x,y)*“^.  (37) 


Integrating  the  resulting  equations  with  respect  to  y  yields  a  wave  equation 
for  pi(x)  in  terms  of  dTm/dx  and  material  properties  and  geometiy.  and  an  equation  for 
energv-  flux  H2  (which  includes  both  heat  and  work)  along  x  tn  terms  of  pi(x),  Tm(x), 
material  properties,  and  geometry.  The  wave  equation  is 

{l+(Y-l)fK/(l+es))pi  +  (pmc2/ci)2)d/dxll(l-fv)./pmldpi/dx) 


-  (pc2/(o2)l(f^fv)/((l-G)(l+es))l(dTm/dx)(dpi/dx)  =  0,  (38) 

where 

fv  =  tanhl(l+l)yo/S\'l/{(l+l)>'o/6v).  (39) 

fK  =  tanhlll+l)yo/6Kl/l(l+l)yo/6K-}.  (40) 

ts  =  pmCp6((tanhl(l+l)yo/6Kl/psCs6stanhl(l+l)L/6s),  (41) 

6v  =  \(2v/(d).  (42) 

5k-^(2k/w),  (43) 
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and 


8s  =  \(2ks/o)). 


(44) 


The  term  6s  is  the  ihermal  penetration  depth  In  the  solid,  and  Is  analogous  to  the 

thermal  penetration  depth  in  the  fluid  discussed  previously,  whUe  6v  is  the  viscous 
penetration  depth  of  the  fluid.  The  viscous  penetration  depth  Is  roughly  the  distance 
shear  diffuses  In  a  time  l/w.  The  equation  for  energy  flux  is 

H2  =  (nyo/2a)pm)Im[(dpi*/dx)pi(l-fv*-Tin  P(fic-fv*)/((l-Kj)ll+Es))ll 

+  (nyocp/2o)3pm(  1  -o))(dTni/dx)(dpi  /dx)(dpr/dx) 

X  ImIfv%|(fvr-fv’)(l+es(fv/fK)/((l+a)(l+es)l]l 

-  ri(v'oK+LKsl(dTm/dx).  (45) 

In  order  to  make  contact  with  the  single  plate  results,  we  make  three 
assumptions  designed  to  simplify'  the  equations.  These  assumptions  are  zero  viscosity, 
the  "boundary-layer"  approximation,  and  the  "short-stack"  approximation.  (The 
assumption  that  the  fluid  viscosity  Is  zero  Is  not  a  good  approximation  In  most  realistic 
situations,  requiring  that  the  viscous  terms  b«.  ut  back  Into  the  equations  at  a  later 

lime.)  With  the  fluid  viscosity  set  equal  to  zero  the  terms  a  and  fv  both  vanish. 

The  ' boundary-layer '  approximation  consists  of  the  assumptions  that 

yo»8K  and  that  L>>6s.  The  goal  of  these  assumptions  Is  to  set  the  hj-pcrbollc  tangents 
equal  to  unity.  This  Is  reasonable  for  cases  where  the  value  of  the  dimension  of  Interest 

Is  at  least  twice  die  value  of  the  associated  boundary  layer  (l.e.  yo  i  25k  llSwlft.  1988) 

The  ".shon-stack"  approximation  assumes  that  the  slack  Is  short  enough 
that  It  does  not  appreciably  disturb  the  standing  wave,  so  that  we  can  lake 

pi  =  Po  co.slkx)  =  pi®  (46) 
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and 


ul  =  -l(l+L/yo)(Po/pmc)sln(kx)  =  -lui®.  (47) 

These  are  the  same  expressions  used  In  the  single  plate  calculations  except  for  the  teitn 
1+L/yo  which  Is  required  for  continuity  of  volume  velocity.  It  Is  further  assumed  that 

the  stack  is  short  enough  that  pi®  and  ui®  may  be  regarded  as  Independent  of  x  within 
the  stack,  and  that  AT=Th-Tc  «Tni.  so  that  thermophysical  properties  can  be  assumed 
Independent  of  Tm  and  hence  of  x. 

With  these  assumptions,  and  using  the  equation 

dpl/dx= -copmul®,  (48) 

(Swift,  1988,  p.  1158)  we  obtain 

H2  =  -(l/4)n5K(TmPpi®ui®/(l+e9)lir-l) 

-  n(yoK+LKs)(dTm/dx)  (49) 

for  energy  flux,  and 

pi  +  (pmc2/u)2)d/dx{(l/pm)dpi/dx] 

=  l(7-l)5Kpi/|(l+l)(l+Es)yo)l(r  -  1)  (50) 

for  the  wave  equation. 

Tlie  acoustic  power  generated  (or  absorbed)  In  the  stack  may  be  determined 
by  multiplying  the  change  In  average  acoustic  intensity  between  the  two  ends  of  the 
stack  by  the  fluid  cross  section  at  the  stack  exit.  The  resulting  equation  is 

W2  =  nyoI<piUl>rlght  -  <piui>lefi] 

=  nyoAxd(<piui>)/dx.  (51) 

Manipulating  Eq.  51.  using  Eqs.  48  and  50.  we  find  that 


W2  =  (l/4)n5KAx|(>-l)a)/(pmc2(l-t-cs)}I(piS)2(r-l) 


(52) 


We  pause  now  to  make  contact  with  the  simple  results  of  the  single  plate 
calculations.  We  begin  by  examining  the  energy  flux  equation.  The  first  term  on  the 

right  in  Ekj.  49  dlflers  from  the  right  side  of  Ek}.  19  by  a  factor  l/(l+es)  which  accounts 
for  the  finite  heat  capacity  of  the  plates  and  modifies  the  boundary  condition  for 
temperature  at  the  fluid-solid  Interface.  The  second  term  on  the  right  side  of  Eq.  49 
accounts  for  ordinary  heat  conduction  in  the  plates  and  the  fluid.  Upon  examination  of 
the  equations  for  work  flux  (Eqs.  25  and  52)  we  find  that  the  difference  Is  again  the 

factor  l/(l+es)  described  above. 

We  now  Include  v-iscoslty  and  repeat  the  above  calculations,  but  now  use  the 
relationship 


dpl/dx  =  -apinui^/d-fv)  (53) 

(Swift,  1988,  p.  1 159).  Notice  that  dpi  /dx  Is  now  complex  due  to  the  presence  of  viscosity 
in  the  fluid.  (Note  that  Elqs.  46  and  47  may  still  be  used  for  this  case  as  long  as  we  now  let 

ui®  represent  the  mean  fluid  velocity  in  the  x  direction.  Due  to  the  presence  of  viscosity, 
the  X  component  of  the  fluid  velocity  Is  a  function  of  y.  going  from  zero  at  the  surface  of 
the  plates  to  some  maximum  value  at  a  distance  on  the  order  of  a  viscous  penetration 
depth  away  from  the  plates.)  We  now  obtain 

H2  =  - 1  /4n  SxITmPp l®u  i®/((  1  +c)(  1  -t-es)!  1  •6v/yo+6v2/2yo2))] 

X  ini+vo+o+aesl/d+vo)  -  (l+vo-6v/yo)l 

-  n(yoK+LKs)(dTm/dx)  (54) 

for  energy  flux,  and 
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pi  +  (pmc2/a)2)d/dx{((l-fv)/pmldpi/dx) 

=  (7-l)6KPl/((l+l)(l+ea)yo}ir/((lWo)(l-fv))  -  1)  (55) 

for  the  wave  equation. 

The  work  flux  is  calculated  In  the  Scime  manner  as  before  and  the  result  Is 
W2  =  (l/4)n5KAx((Y-l)a>(pi®)2/{pnic2(l+es))l 
X  |r/((l+vo)(l-5v/yo-hSv2/2yo2))  -  1] 

-  (l/‘.)n&vAxla)pm(uiS)2/(i.5v/yo+6v2/2yo2).  (56) 

It  is  clear  that  the  Inclusion  of  viscosity  in  the  calculations  adds  significant  complexity 
to  the  results,  and  the  intuitive  understanding  of  the  components  of  the  equations  is 
lost. 
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m.  DESCRIPTION  O  '  APPARATUS 


The  apparatus  which  we  are  modeling  In  this  thesis  is  a  thermoacoustic  prime 
mover,  enclosed  in  a  rigid  tube  resonator,  driving  an  acoustlc-to-electrlc  power 
transducer.  Figure  8  is  a  general  schematic  of  the  device.  Specific  values  of  the 
dimensions  of  the  prime  mover  are  contained  Ln  Table  1 . 

The  transducer  proposed  for  use  In  this  application  is  a  linear  alternator.  Design 
requirements  submitted  for  bid  by  Dr.  Greg  Swift  of  L.A.N.L.  Include  1  kW  of  power 


Device 

Characteristic 

Hot  heat 
exchanger 

Thermo - 

acoustic 

stack 

Cold  heat 
exchanger 

Inlet  position 
relative  to  rigid 
end  (m) 

0.305 

0.356 

0.584 

Ebdt  position 
relative  to  rigid 
end  (m) 

0.356 

0.584 

0.635 

Plate  thickness  (m) 

4.57e-4 

1  27e-4 

2.54e-4 

Fluid  gap  (m) 

9.40e-4 

7.94e-4 

7.94e-4 

Material 

Ni 

304  SS 

Cu 

K(W/mK) 

89.9 

16.3 

350.0 

Cp(J/KgK) 

443.8 

451.9 

384.0 

Density 

(Kg/m^3) 

8890 

8027 

8900 

Table  1:  Summary  of  geometry  of  thermoacoustic  stack  assembly  am 
important  thermodynamic  parameters. 

delivered  to  an  electric  load  resistance  of  3  to  100  ohms  with  at  least  60%  efficiency 
when  driven  by  a  4  bar  peak-to-peak  pressure  amplitude  at  120  Hz  in  helium  at  a  mean 
pressure  of  20  bar.  The  details  of  the  proposals  received  in  response  to  the  request  for 
quotation  contain  proprietary  Information,  The  transducer  model  used  In  this  thesis 
was  designed  with  the  assistance  of  Nick  Vitale  of  MTI.  Although  It  does  not  represent 
the  actual  design  proposal  submitted  by  MTI,  It  should  provide  reasonable 
approximations  for  the  performance  of  the  prime  mover-transducer  system.  The 
transducer  may  be  modeled  by  the  equivalent  circuit  Indicated  In  Fig.  9,  The  moving 

mass  of  the  transducer  is  0.766  Kg,  the  sprtng  constant  is  1,43E5  Kg/s^,  the  magnetic 
field  strength  Is  0.97  Tesla,  and  the  coll  has  a  total  effective  length  of  32.6  m.  The 
mechanical  resistance  of  the  transducer  Is  2.0  Kg/s.  Although  the  actual  transducer  will 
have  some  physical  limitation  to  its  displacement,  we  have  not  Imposed  any  such 
limitation  on  our  model. 
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Figure  9:  Equivalent  circuit  model  of  linear  alternator  acoustlc-to-electric  power 
transducer. 


IV.  NUMERICAL  MODEL 


In  order  to  formulate  a  program  capable  of  modeling  the  peiformance  of  our 
thermoacoustic  generator,  we  must  Identify  the  parameters  which  may  be  allowed  to 
vary  during  the  operation  of  the  device,  and  then  develop  a  system  of  equations  that 
will  relate  these  variables.  The  geometry  of  the  device  will  be  fixed.  It  Is  assumed  that 
the  thermodynamic  properties  of  the  material  used  in  the  construction  of  the  device 
will  not  vary  significantly  over  the  range  of  pressures  and  temperatures  at  which  the 
device  will  operate.  It  is  further  assumed  that  the  temperature  of  the  cold  heat 
exchanger  will  be  constant.  The  working  fluid  used  will  be  a  mixture  of  helium  and 
argon  gases.  The  mixture  of  these  gases  will  be  determined  so  as  to  achieve  a  120  Hz 
operating  frequency  at  a  median  load  resistance  of  13  ohms.  The  parameters  that  will 
be  allowed  to  vary  in  the  model  are  the  acoustic  pressure  ampbtude.  angular  frequency, 
and  the  temperature  of  the  hot  heat  exchanger,  for  a  given  mput  heat  flux  and  load 
resistance.  Variation  of  the  hot  end  temperature  will  result  in  variation  of  the 
themiodynamlc  properties  of  the  working  fluid  and  of  the  sound  speed  In  the  fluid. 

Before  continuing,  one  additional  assumption  will  be  made  concerning  the 
performance  of  the  system;  specifically,  that  no  attenuation  occurs  in  the  resonator. 
All  acoustic  power  generated  in  the  stack  must  therefore  be  removed  at  the  transducer. 
This  assumption  is  made  to  simplify  the  calculations,  but  is  not  truely  valid  since  a 
significant  fraction  of  the  acoustic  power  would  be  attenuated  due  to  thermal  and 
viscous  effects  at  the  walls  of  the  resonator. 

We  now  begin  to  develop  the  equations  which  will  determine  the  performance  of 
the  system  using  the  short  stack  and  boundaiy  layer  approximations  discussed 
previously.  Starting  at  the  rigid  end  of  the  resonator,  and  assuming  a  standing  wave,  we 
know  that  the  pressure  amplitude  must  be  at  a  meixlmum  and  the  fluid  particle  velocity 
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must  be  zero.  We  can  therefore  describe  the  pressure  and  fluid  particle  velocity  in  the 
region  between  the  rigid  end  and  the  stack  assembly  using  Eqs.  6  and  8.  which  are 
repeated  below  for  ease  of  reference. 


p(x)  =  Po  cosdoc) 

(6) 

u(x)  =  -l(Po/pmC)  sln(kx) 

(8) 

The  short  stack  approximation  allows  us  to  assume  that  the  pressure  amplitude  in  the 
stack  Is  constant.  Using  the  mean  stack  position  to  determine  the  value  of  this  term  we 
find 


p(xm)  =  pi®  =  Po  cos(kmxm) 

(57) 

where 

km  =  co/cm. 

(58) 

and  Cm  is  the  mean  sound  speed  in  the  system.  The  mean  fluid  particle  velocity  at  the 
mean  stack  position  is  given  by 

ufxrrO  =  -lui®  =  -lg(Po/pmc)  sln(kmXm),  (59) 

where  g  =  1+L/yo  is  the  geometry  factor  necessary  to  maintain  continuity  of  volume 

velocity.  It  must  be  noted  that  g  has  different  values  for  the  slack  and  the  two  heat 
exchangers. 

The  change  In  acoustic  pressure  due  to  the  stack  and  heat  exchangers  Is  obtained 


from  Eq.  53.  The  resulting  equations  are 

dpht  =  -PoWdxhtght  sln(kmXm)/lCm(l-fv.ht)l.  (60) 

dpst  =  -Poudxstgst  sln(kmXni)/lcm(Tfv,st)l.  (61) 

and  dped  =  ■PowdXfdgrd  sln(kmXm)/{cm(Tfv,cd)}.  (62) 

where  fv  =  1 /((I +i)yo/8v).  (63) 


Note  the  subscripts  indicating  the  parameters  which  have  different  values  depending 
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on  the  location  In  the  stack  assembly.  The  acoustic  pressure  at  the  exit  of  the  stack 
assembly  is  then  the  sum  of  the  pressure  changes  in  the  stack  and  the  pressure  at  the 
stack  entrance: 


p(x4)  =  p(xi)+dpht+dpst-Kipcd 

(64) 

where 

p(xi)  =PoCos(khtXl). 

(65) 

The  change  in  mean  fluid  particle  velocity  Is  obtained  from  the  wave  equation 
(Eq.  55)  by  substituting  for  dpi/dx  from  Eq.  53  and  solving  for  du.  The  resulting 
equation  is 

du  1 =l[o3p  1  dx/  (pmc'-^)l 

X  (((Y-l)SK/((l+l)(l+Es)yo))(r/((l+Vo)(l-fv])-ll-ll.  (66) 

The  changes  in  mean  fluid  particle  velocity  In  the  heat  exchangers  and  the  stack  are 
then 

duht=il«Podxhl  cos(kmXm)/(pmCm2)) 

X  l((7-l)5K/{(l+l)(l+es,ht)yo.ht))ir/((li-va)(l-fv.hO)-ll-ll.  (67) 

dust=llwPodxst  cos(kmXm)/(pmCm2ll 

X  (((Y-l)SK/((l+l)(l+Es,sl)yo.st))ir/((l+Vo)(l-fv.st))-ll-ll.  (68) 

and 

ducd=l((oPodxcd  cos(kmXtn)/(pmCm'‘i! 

X  1((7- 1  )5k/((  1  +1)(  1 4E5,cd)yo,cd))ir/((  1  +vo)(  1  -fv.cd))  - 1 1- 1 ).  (69) 

The  fluid  particle  velocity  at  the  stack  exit  Is 

u(x4)  =  u(xi)+duht/ght+dusl/gst+ducd/gnd  (70) 

where  u(xi)  =  -i(Po/rmCht)  sln(khixi).  (71) 
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The  geometry  factors  are  required  for  conttnult>'  of  volume  velocity  as  before. 

Assume  now  that  the  acoustic  wave  in  the  region  between  the  stack  assembly  and 
the  transducer  has  the  form  of  left  and  right  traveling  plane  waves  described  by 


plx.t)  ^  (72) 

where  A  and  B  eire  complex  coefficients  given  by 

A  =  a  +  la  (73) 

and  B  =  b  +  l3.  (74) 

The  associated  fluid  particle  velocity  is  found  from  Euler's  equation  to  be 

u(x,t)  -  tA/pnjCcd)e*(“^‘^  -  (B/pmCcd)e‘l“^‘*^^.  (75) 


Dropping  the  time  dependence  of  Eqs.  72  and  75.  and  expanding  In  terms  of  sines  and 
cosines  yields  equations  for  acoustic  pressure  and  fluid  particle  velocity  In  the 
resonator  between  the  slack  assembly  and  the  transducer  of  the  form 

p(x)  =  (a+b)cos(kcdx)+(a-p)sln(kcdx) 

+  ll(a+(l)cos(kcdx)-(a-b)sln(kcdx)l  (76) 

and 

u(x)  =  (l/pniCcd)((a-b)cos(kcdx)+(a+P)sln(kcdx) 

+  ll(a-  P)cos(kcdx)-(a+b)sln(kcdx)l)  (77) 

One  set  of  boundary  conditions  for  Eqs.  76  and  77  are  given  by  our  expressions  for 
pressure  and  velocity  at  the  stack  assembly  exit.  This  requires  expansion  of  the 
expressions  into  real  and  Imaginary  parts  so  that  we  may  form  the  equalities 

(a+b)cos(krdx)+(a-p)sln(kcdx)  =  Re{p(xi)+dphi+dpst+dpcd),  (78) 

(a+p)cos(kcdx)-(a-b)sln(kcdx)  =  lm(p(xi)+dphl+dpsl+dpcd).  (79) 

(1  /pmccd)i(a-b)cos(kcdx)+(a+P)sln(kcdx)l 
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=  Re(u(xi)+duht/ght+dust/gst+ducd/gcd}. 


(80) 


and  ( 1  /pmCcd)((a-p)cos(kcdx)-(a+b)sin{kcdx)) 

=  lm(u(xi)+duhl/ght+dust/gst+ducd/gccl}.  (81) 

The  boundary  condition  at  the  transducer  required  for  resonance  Is  that  the 
reactance  of  the  combined  system  must  vanish: 

Im{Zr  +  Zm)  =  0.  (82) 

The  mechanical  impedance  of  the  acoustic  wave  at  the  transducer  position  is  given  by 

Zr=(7td2/4)lp(x0/u(xt)l.  (83) 

The  mechanical  impedance  of  the  transducer  is  found  by  applying  Newton  s  second  law 
to  the  transducer.  The  force  acting  on  the  transducer  by  the  acoustic  pressure  Is 


Ft  =  (Ttd2/4)p(xt)  (84a) 

Applying  Newton's  second  law  to  the  transducer  gives 

Ft  =  mia  +  +  Flmu +BLcl.  (84b) 

This  equation  may  be  ref  ist  by  recognizing  that 

a  =  iiou,  185) 

;  =  -lu/G).  (86) 

and  1  =  BLcu/Ze.  (87) 

where  Ze  =  1 1  /Rl  + 1  /{R2+Rload+l(o>L- 1  /wC)})'  ^ .  (88) 

The  result  is  that  the  mechanical  Impedance  of  the  transducer  may  be  written  as 

Zm  =  Ft/u  =  Rm  +(BLc)2/Zc  +l(tijmi-s/o)).  (89) 


We  may  find  two  more  equations,  necessary  to  fully  constrain  the  7  variables  of 
Interest,  by  recognizing  that  the  conditions  in  the  system  are  directly  effected  by  the 
amount  of  heat  flux  in  the  stack,  and  that  the  power  propagating  from  the  stack 
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assembly  must  be  removed  at  the  transducer  during  steady  state  operation.  Our  sixth 
equation  Is  therefore  Elq.  54.  The  seventh  equation  is 


riac  -  rim  -  ric  =  0  (90) 

where  Oac  =  (JKl^/8)Re(pu*}.  (91) 

nm  =  (l/2)Rm{uu*}.  (92) 

and  ne  =  (l/2)Re(I*V).  (93) 

The  voltage  produced  by  the  transducer  Is 

V  =  BLcu.  (94) 

Before  describing  the  method  used  to  Arid  solutions  to  these  equations,  It  Is  useful 
to  summarize  them  and  the  variables  of  Interest.  The  seven  equations  are 

(1)  (a+b/cos(k£dx)+(a-P)slri(kcdx)  =  Re|p(xi)+dphH  dpsi+dpcd), 

(2)  (a-(-P)cos(kcdx)-(a-b)sin(kcdx)  =  lm!p(xi)+dphi+dpst+dpcdl. 

(3)  (l/pmCcd){(a-b)cos(kcdx)+(u+P)sln(kcdx)) 

=  Re(u(xi)+duhl/ght-Klust/giil+ducd/gcd). 


(4)  (l/pniCcd){(a-p)cos(kcdx)-(a+b)sln(ki:dx)) 

=  lm(u(xi)+duhi/ghi+dust/gsi+ducd/gcd). 


(5)  lm(Zr  +  Zm)  =  0  . 

(6)  H2  =  -l/4n5K|TmPpi®uiVi(l+o)(l+£s)(l-6v/yo+6v2/2yo2))l 

X  |ri(l+\a+o+a£s)/(( l+vo))  -  (l+Vo-6v/yo)l 
-  n(yoK+LKs)(dTm/dx), 
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and 

(7)  riac  *  rim  -  He  =  0. 

The  variables  for  which  we  are  to  solve  these  equations  simultaneously  are  the  five 
pressure  coefficients  Po.  a,  a,  b,  and  p,  the  angular  frequency  ai,  and  the  temperature  of 
the  hot  heat  exchanger  Tht.  In  order  to  accomplish  this  task  effectively,  it  Is  necessary 
to  expand  the  above  equations  In  terms  of  the  parameters  of  Interest  and  to  solve 
analytically  for  the  real  and  Imaginary  parts  of  the  equations  Indicated,  This  was 
accomplished  using  an  algebraic  manipulation  software  system  called  Derive;  a 
Mathematical  Assistant  for  Personal  Computers.  The  details  and  results  of  the 
calculations  are  to  long  to  present  here. 

Before  we  begin  to  solve  these  equations  we  must  know  the  values  of  all  the 
parameters  that  are  not  variables.  Tlie  parameters  that  describe  the  geometry  of  the 
system  and  the  characteristics  of  the  load  will  be  fixed  at  the  beginning  of  the  process. 
The  parameters  that  describe  the  working  fluid  are  dependent  on  the  gas  mix  used  and 
on  the  temperature  at  the  hot  heat  exchanger.  This  temperature  effects  the  values  of  the 
viscosity,  thermal  conductivity,  and  Prandtl  number  of  the  gas  mix  as  well  as  the, 
sound  speed  in  the  s>’stem.  which  in  turn  effects  the  gas  mix  that  will  be  chosen  in  order 
to  have  resonance  at  the  desired  frequency.  It  is  necessary  therefore,  to  begin  the 
calculations  by  guessing  a  value  of  the  hot  end  temperature  and.  based  on  this  guess, 
calculate  the  required  gas  p.irameters.  We  must  then  solve  the  equations  to  find  a  new 
hot  end  temperature,  and  then  adjust  the  values  of  the  gas  parameters  In  an  Iterative 
cycle,  until  a  self-consistent  solution  Is  obtained. 

We  begin  the  process  by  "guessing  "  a  value  for  Tht-  Using  this  value  we  calculate  the 
mean  temperature  gradient  VTm.  the  mean  temperature  In  the  slack  Tm,  the  average 
system  temperature  Tavg.  and  the  critical  temperature  gradient  VTcrlt  using  the 
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relations 


VTm  =  frht-Tcd)/dxst. 

(95) 

Tm  =  (Tht+Tcci)/2. 

(96) 

Tavg  =  (xaTht  +dxstTni  +(x5-X3)Tcd)/x5. 

(97) 

and 

VTcrlt  =  I(Y*l)/Tmp]{Tmkm)  col(knixin), 

(17) 

where 

k  =  OJ/cm  =  it/xt. 

(98) 

For  this  first  approximation  it  is  assumed  thatTmP  =  )  and  y=1.67.  Using  the  result  of 
the  above  calculations  we  then  determine  the  normall2ed  temperature  gradient 

r=VTm/VTcr1t. 

With  the  temperature  parameters  calculated  above  we  now  determine  the  first 
estimate  of  the  gas  mix  required  to  obtain  resonance  at  the  desired  frequency  fo.  This  is 
accomplished  by  using  the  expression  fc.  sound  speed  in  a  perfect  gas 

c  *=  Vlyr'D  (99) 

and  estimating  the  required  sound  speed  as  c  =  foX  =  2foxt.  The  gas  constant  is  given  by 

r  =  R/(FheWhe  +  FarWar)  1 100) 

where  R  is  the  universal  gas  constant.  F  is  the  mole  fraction  of  the  gas  indicated  by 
subscript  that  is  present,  and  W  is  the  molar  weight  of  the  gas  Indicated  by  subscript. 
Since  only  the  two  gases  are  present  in  the  system,  we  may  write 

Far=l-Fhe.  (101) 

Substituting  Eqs.  100  and  10)  Into  Eq.  99.  and  using  c  =  2foxt.  we  solve  for  the  mole 
fraction  of  helium  required: 

Fhe  =  1  /(Whe-War)lYRTavg/(2foXt)2  -  War).  ( 102) 
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From  these  results  we  c£ilculate  the  molar  weight  of  the  binary  gas  mixture  for  future 
use: 

Wmix  =  FhcWhe  +  FarWar.  ( 103) 

Having  established  an  approximation  for  the  required  gas  mix,  we  next  determine 
the  Isobaiic  and  isochorlc  specific  heat  capacities  of  the  gases  at  Tavg.  This  is  done  by 
two  subroutines  called  'HELIUM'  and  ‘ARGON’.  These  subroutines  are  modlflcatloirs  of 
programs,  written  by  Robert  D.  Mccarty  at  the  National  Bureau  of  Standards,  to 
calculate  the  thermodynamic  properties  of  these  and  other  gases,  using  a  form  of  the 
vliial  equation  of  state.  The  modifications  to  the  programs  consisted  primarily  of 
eliminating  functions  not  required  for  this  application,  such  as  output  to  files  and 
calculation  of  additional  thermodynamic  properties  (Internal  energy,  entropy, 
enthalpy,  etc..),  and  modifying  the  input  requirements  so  that  the  subroutines  read 
temperature  and  pressure  (the  Input  parameters  required)  from  the  main  program 
instead  of  from  the  screen  In  an  Interactive  mode.  In  addition,  the  upper  temperature 
limit  In  'ARGON'  was  raised  from  400  K  to  10(X)  K.  Comparison  of  data  calculated  in  this 
extended  temperature  region  was  found  to  be  in  agreement  with  tabulated  data 
contained  In  available  tables  of  thermodynamic  properties  for  argon.  Using  the  output 
from  the  subroutines,  the  Isobarlc  and  Isochorlc  specific  heat  capacities  of  the  binary 
gas  mix  are  deteimlned: 

Cp.inlx  =  (cp.heFheWhe  +  Cp.arFarWar)/Wmlx  (104) 

and  Cv.mlx  =  (cv.heFheWhe  +  Cv,arFarWar)/Winlx.  (105) 

From  these,  we  then  calculate  a  new  value  of  the  ratio  of  heat  capacities,  y  = 
cp.mlx/cv.mlx.  If  Y  varies  significantly  from  the  Initial  value  of  1.67,  the  gas  mix  Is 
calculated  again,  using  this  new  value,  until  a  self  consistent  value  is  achieved. 

Now  that  a  reasonable  gas  mix  has  been  established,  and  a  value  of  the  ratio  of 
heal  capacities  determined,  we  can  calculate  the  mean  gas  density  and  acoustic  sound 
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speed  In  the  resonator,  assuming  Ideal  gas  behavior.  The  meem  gas  density  Is  found 
using  the  Ideal  gas  relation 

Pm  =  pmWmtx/RTavg.  (106) 

The  acoustic  sound  sp>eed  in  the  resonator  Is  found  using  Elq.  99. 

Of  more  interest  are  the  gas  properties  in  the  stack  assembly,  where  the 
thermoacoustic  effects  are  occurring.  To  determine  these,  we  again  use  the  subroutines 

HELIUM'  and  'ARGON',  but  with  Tm  as  the  Input  temperature.  We  obtain  the  specific 
heat  capacities  of  the  two  gases,  calculate  the  specific  heat  capacities  of  the  binary 
mixture,  and  determine  y  and  the  acoustic  sound  speed  for  the  stack  assembly  as 
described  above  for  the  resonator.  To  determine  the  dynamic  viscosity,  thermal 
conductivity,  and  Prandtl  number  of  the  binary  gas  mixture,  we  call  subroutine 
PRANTL'.  PRANTL  Is  a  modified  version  of  a  program  written  by  a  student  at  the  Naval 
Postgraduate  School  (Susalla.  1988)  to  calculate  these  same  parameters  for 
thermoacoustic  refrigerator  design  calculations.  The  original  program  calculated  the 
dynamic  viscosity,  thermal  conductivity,  and  Prandtl  number  for  helium-argon  and 
hellum-.xenon  gas  mixes  at  only  two  temperatures.  The  program  was  modified  to  allow 
for  calculations  over  a  range  of  temperatures,  but  only  for  helium-argon  mixtures.  The 
input  parameters  required  for  the  subroutine  are  the  temperature  of  interest  and  the 
Isobarlc  specific  heat  capacities  of  the  two  gases. 

Having  now  obtained  the  necessary  parameters  based  on  our  first  estimate  of  the 
hot  end  temperature,  we  arc  nearly  ready  to  solve  our  seven  equations.  The  seven 
equations  which  we  must  solve  simultaneously  are  non  linear  In  the  variables  of 
interest,  and  can  not  easily  be  solved  analytically.  We  chose  Instead,  to  use  a  readily 
available  Fortran  Library  subroutine  to  find  the  solution  for  us.  The  subroutine  used  Is 
the  IMSL  Library  subroutine  NEQNF.  NEQNF  Is  a  smgle  precision  subroutine  which 
uses  The  Levenberg-Marquardt  algorithm  and  a  flnlle-dirference  approximation  to  the 
Jacobian  to  solve  a  system  of  nonlinear  equations. 
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In  order  to  find  the  solution,  the  subroutine  requires  Initial  "guesses"  for  the 
values  of  the  variables.  We  have  established  a  "guess  "  for  the  hot  heat  exchanger 
temperature  in  order  to  estimate  the  gas  properties  above.  The  ’"guess’"  for  the  angular 

frequency  is  obviously  o)  =  27ifo.  What  is  lacking  arc  Initial  values  of  the  pressure 
coefficients. 

In  previous  work  with  the  5  Inch  engine  at  L.A.N.L.,  operating  it  in  a  simple  rigid- 
rigid  resonator,  acoustic  pressures  on  the  order  of  10%  of  the  mean  pressure  have  been 

observed.  We  will  take  this  value  for  Po.  Values  for  a.  a,  b,  and  p  may  be  obtained  by 
analytically  solving  equations  similar  to  Eqs.  78-81  derived  from  the  Invlscld 
relationships  for  the  short  engine.  The  resulting  equations  are  too  long  to  be  repeated 
here.  Initial  ’’guess'"  values  for  a,  a,  b.  and  p  are  calculated  in  subroutine  "FIRST. 

Having  determined  all  of  the  required  input  parameters,  we  now  call  the  library 
subroutine  to  solve  our  system  of  equations  for  the  seven  "unknown"  variables.  If  the 
resulting  new  hot  end  temperature  differs  slgnlffcantly  from  the  initial  "guess"  It  Is 
necessary  to  redo  the  calculations  of  the  temperature  dependent  parameters  and  gas 
parameters  (excluding  the  gas  mix)  and  solve  again.  This  process  is  repeated  until  the 
values  of  the  hot  heat  exchanger  temperature  do  not  change  significantly  from  one 
Iteration  to  the  next. 

Once  a  consistent  value  of  Tht  Is  determined,  we  compare  the  resulting  angular 

frequency  to  our  desired  operating  frequency.  If  the  calculated  frequency  Is  higher  than 
the  design  point,  we  increase  the  argon  fraction,  thereby  reducing  the  acoustic  sound 
speed,  and  repeat  the  iteration  process,  Conversly,  If  the  frequency  is  too  low,  the  argon 
fraction  is  reduced,  and  the  Iteration  process  repeated.  These  corrections  are  repeated 
until  a  gas  mix  is  found  which  puts  use  near  the  desired  operating  frequency. 

Finally  we  compare  the  electric  power  generated  by  the  transducer  (not  that 
delivered  to  the  load)  with  a  desired  maximum  value.  If  the  power  generated  differs 


38 


from  the  design  value.  H2  is  adjusted  and  the  entire  Iteration  process  Is  repeated. 

The  ultimate  result  of  the  above  process  is  a  gas  mix  that  will  result  in  operation 
of  the  system  at  the  desired  frequency,  and  provide  the  desired  power  output,  for  a  given 
load  reslstEince  and  a  particular  heat  flux  input. 

It  should  be  noted  here  that  the  above  iterative  processes  are  actually  carried  out 
In  two  stages.  The  first  stage,  accomplished  In  subroutine  'INITIAL',  provides  a  rough 
estimate  of  the  required  values.  Subroutine  'REFINE'  takes  these  rough  estimates  and 
refines  them  to  the  final  values. 

Once  the  gas  mix  is  determined,  it  Is  fixed  and  the  system  of  equations  may  be 
solved  repeatedly  for  a  variety  of  Input  heat  fluxes  and  load  resistances  to  determine 
the  response  of  the  system  about  the  design  point.  The  instantaneous  power  delivered  to 
the  load  by  the  transducer  is  given  by 

rii=VI.  (107) 

The  current  flowing  through  the  load  resistor  is  given  by 

I  =  (V/Z)  sln(ut+4»)  (108) 

(Ohanlan.  1985)  where  V  Is  the  voltage  produced  by  the  transducer  given  by  Eq.  137.  ^  Is 
the  phase  difference  between  the  voltage  and  the  current  given  by 

<t)  =  tan-ll(l/(i)C-coL)/(R2  +  Rload)).  (109) 

and  Z  Is  the  Impedance  of  the  branch  containing  the  load  resistor  given  by  the  relation 
Z  =  R2  +  Rjoad-l(0JL-l/uC).  (110) 

After  appropriate  substitution  and  reduction  we  find  that  the  time  average  power 
delivered  to  the  load  Is 

nioad  =  l/2(BLc)'^Rjoad{u'u/Z*Z)  cos((t)).  (Ill) 

The  term  cosW  Is  very  nearly  unity  for  the  values  of  L.  C.  and  cu  of  interest,  and  Is 
neglected  in  the  calculations  of  Hload  in  the  program.  We  may  now  ceilculate  the  total 
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system  efficiency  as 


fi  =  nioad/H2. 


(112) 
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V.  RESULTS 


In  this  chapter  we  will  present  the  results  of  calculations  made  for  various  load 
resistances,  input  heat  fluxes,  and  transducer  positions.  The  trends  of  some  Important 
thermoacousllc  variables  are  presented  graphically.  Since  no  experimental  device 
comparable  to  the  one  modeled  has  been  built,  no  comparison  of  theoretical  and 
experimental  data  is  possible.  The  complexity  of  the  governing  relationships  makes 
Intuitive  prediction  of  the  trends  dllTlcult. 


A.  VARIATION  IN  LOAD  RESISTANCE  AND  INPUT  HEAT  FLUX 

In  order  to  Investigate  the  effect  of  varying  the  load  resistance  and  input  heat 
flux,  we  first  fix  the  transducer  position  at  4.0  m.  We  then  calculate  the  gas  mix 
required  for  resonance  at  120  Hz.  with  a  nominal  load  resistance  of  13.0  O.  The  result 
of  these  calculations  is  a  gas  mix  of  97%  helium  and  3%  argon. 

Table  2  contains  some  representative  values  of  important  parameters  calculated 


Load 

Resistance 

(fl) 

Angular 

Frequency 

(rad/s) 

Hot 

End 

Temp 

(K) 

Useful 

Power 

Generated 

(W) 

Overall 

Efficiency 

(%) 

Acoustic 

Pressure 

Amplitude 

(Bar) 

B 

12 

9 

766 

1100 

589 

11.8 

1.26 

10 

764 

1060 

577 

11.5 

1.29 

11 

762 

1030 

565 

11.3 

1.31 

12 

760 

1000 

552 

11 

1.33 

13 

758 

974 

528 

10.8 

1.34 

14 

757 

951 

527 

10.5 

1.37 

15 

755 

931 

514 

10.3 

1.4 

16 

754 

913 

502 

10 

1.41 

17 

753 

897 

491 

9.8 

1.44 

Table  2:  Values  of  critical  parameters  as  a  function  of  load  resistance  for  a  97% 
hebum  gas  mix.  The  transducer  is  at  4.0  m  and  the  input  heat  flux  Is  constant  at  5 
kW. 
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for  an  Input  heat  flux  of  5  kW.  Figures  10  through  14  Indicate  graphically  the  general 
trends  of  these  parameters  as  the  load  is  increased.  Figure  15  shows  the  variation  in 
acoustic  pressure  amplitude  as  a  function  of  input  heai  flux  for  a  number  of  load 
resistances.  The  output  of  the  program  indicated  no  variation  in  frequency,  hot  end 
temperature,  or  efficiency  as  the  liiput  heat  flux  was  changed  ^  '  a  given  load  resistance 
so  no  plots  of  these  relationships  are  presented. 


&  VARIATION  IN  TRANSDUCER  POSITION 

Calculations  were  also  done  for  transducer  positions  of  3.5.  3.0.  and  2.5  meters 
Table  3  and  the  figures  which  follow  Indicate  the  general  trends  in  the  parameters  of 
Interest  as  the  transducer  Is  moved  closer  to  the  Uiermoacoustlc  stack. 

We  clearly  expect  the  argon  percentage  to  Increase  so  that  we  maintain  a  resonant 
system.  It  Is  unclear  what  effect  this  will  have  on  the  overall  performance  of  the 
system.  Changing  the  ratio  of  argon  to  helium  will  certainly  change  the 
thermodynamic  properties  of  the  gas  mix,  but  whether  that  change  will  be  adverse  or 
beneficial  is  not  intuitively  obvious. 

Examining  the  following  figures,  we  see  that  the  normalized  temperature 
gradient  increases,  the  pressure  amplitude  Increases,  and  the  system  efficiency 
decreases  as  the  resonator  is  shortened  for  the  same  stack  geometry. 

The  Increase  In  the  normalized  temperature  gradient  as  the  transducer  Is  moved 
closer  to  the  stack  assembly  may  be  understood  by  noting  that  the  concurrent  change  In 


Xdcr 

Posit 

(m) 

Helium 
Fraction 
(mole  %) 

Thermal 

Penetration 

Depth 

(mm) 

Viscous 

Penetration 

Depth 

(mm) 

Acoustic 

Pressure 

Amplitude 

(Bar) 

Norm. 

Temp. 

Grad. 

System 

Eff. 

(%) 

W1 

97 

1.83 

3.443 

niB 

92 

0.212 

0.152 

1.87 

3.572 

10.4 

Bi 

84 

0.202 

0.138 

1.93 

3.787 

9.66 

70 

0.186 

0.124 

2.03 

4.207 

8.52 

Table  3!  EUect  oi  decreasing  resonator  length  on  parameters  ol  Interest. 
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the  gas  mix  results  In  a  decrease  In  the  sound  speed  and  thermal  penetration  depth  of 
the  gas.  These  decreases  result  In  an  Increase  in  the  temperature  gradient  in  order  to 
produce  the  same  amount  oi  acoustic  power. 
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Figure  13:  Variation  In  system  efficiency  as  a  function  of  load  resistance  at 
constant  heat  flux. 
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Figure  16:  Mole  fraction  of  helium 
13  ohm  load  resistance. 


various  transducer  positions  for 


Figure  17;  Variation  of  penetration  depths  with  helium  fraction. 
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Figure  20:  System  efTiclency  as  a  function  of  transduce-  position. 


» 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


It  Is  not  always  obvious  that  the  trends  indicated  by  the  figures  in  the  proceeding  i 

chapter  are  those  which  we  should  expect.  The  lack  of  a  working  thermoacoustic 
generator  for  comparison  does  not  allow  for  experimental  validation  of  the  results 
obtained.  The  results  are  of  the  correct  order  of  magnitude,  but  the  details  of  the 
relationships  require  additional  Investigation.  It  seems  that  a  reasonable  framework 
has  been  established  for  a  model. 

One  way  In  which  the  program  may  be  Improved  Is  by  Inclusion  of  viscous  and 
thermal  effects  In  the  resonator.  If  this  were  done  then  the  transducer  parameters  could 
be  adjusted  so  that  its  mechanlccil  impedance  would  be  equivalent  to  that  of  a  rigid  end. 

Comparison  of  the  output  of  the  program  to  existing  experimental  data  would  then  be 
possible.  This  rrvay  serve  to  validate  the  output  of  the  program  to  some  degree. 

The  next  requirement  Is  that  a  thermoacoustic  generator  be  built  and  tested.  It 
need  not  be  as  large  as  the  five  inch  engine  at  LjVN.L.  to  measure  general  performance 
trends.  If  the  program  is  modified  as  suggested  above  hen  a  much  smaller  device  in 

* 

which  the  viscous  and  thermal  effects  In  the  resonai  r.re  substantial  may  be  as  easily 
modeled.  A  smaller  device,  working  at  a  mucii  !r  •  „i  pressure  would  not  require  such  an 
expensive  trcinsducer. 


« 
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